Supersonic quantum communication 
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When locally exciting a quantum lattice model, the excitation will propagate through the lattice. The effect is 
responsible for a wealth of non-equilibrium phenomena, and has been exploited to transmit quantum informa- 
tion through spin chains. It is a commonly expressed belief that for local Hamiltonians, any such propagation 
happens at a finite "speed of sound". Indeed, the Lieb-Robinson theorem states that in spin models, all effects 
caused by a perturbation are limited to a causal cone defined by a constant speed, up to exponentially small 
corrections. In this work we show that for translationally invariant bosonic models with nearest-neighbor inter- 
actions, this belief is incorrect: We prove that one can encounter excitations which accelerate under the natural 
dynamics of the lattice and allow for reliable transmission of information faster than any finite speed of sound. 
The effect is only limited by the model's range of validity (eventually by relativity). It also implies that in non- 
equilibrium dynamics of strongly correlated bosonic models far-away regions may become quickly entangled, 
suggesting that their simulation may be much harder than that of spin chains even in the low energy sector. 



Quantum spin chains — or more generally, quantum spin 
models on a lattice — are ubiquitous in condensed matter 
physics and quantum optics. They share the fundamental fea- 
ture that perturbations will propagate through the lattice at 
some characteristic "speed of sound" fljj, |2J] . This effect plays 
an important role for a wealth of non-equilibrium phenom- 
ena in many-body systems, e.g., for the dynamics of relax- 
ation processes towards equilibrium JH0]. In the context of 
quantum information science, it has been noted that excita- 
tions propagating through a spin chain may be used to transmit 
quantum information — thus turning a spin chain into a quan- 
tum channel. Here, the appealing feature is that the transport 
is not facilitated by engineered quantum gates, but rather by 
the natural time evolution of the lattice system lH, eI 0] • 

Because in lattice models only neighboring systems interact 
with each other directly, it is intuitive to assume that the maxi- 
mal propagation speed of excitations (i.e., the speed of sound) 
is finite and given by a value characteristic for each model. 
Indeed, an analogous statement is clearly true for relativistic 
systems, where a perturbation can have no influence outside 
its causal cone. Mathematical physics provides a rigorous jus- 
tification for this observation in the form of Lieb-Robinson 
bounds [1]: in spin lattice systems, perturbations can spread 
only linearly in time, up to exponentially small corrections. 
Recently, an analogous result has been proven to hold for a 
class of bosonic systems 0]. 

Interestingly, familiar as the belief that propagation of ex- 
citations in local models happens with a finite velocity may 
be: it is not quite right. We demonstrate that certain well- 
defined local bosonic models allow excitations to accelerate 
to arbitrarily high velocities. The effect is only limited by the 
range of validity of the model (which must certainly break 
down with the onset of relativistic effects). It occurs even 
for single excitations with bounded energy, traveling along a 
one-dimensional chain of bosons with translationally invari- 
ant nearest-neighbor interactions. From the quantum informa- 
tion perspective, we show that the quantum channel associated 
with this chain has a strictly positive information capacity, 
even after a time sub-linear in the length of the chain. While 
the presented models are non-integrable, we derive the results 



FIG. 1: In quantum spin chains, local excitation will travel no faster 
than some characteristic speed of sound. The phenomenon defines 
a causal cone, outside of which any influence is exponentially sup- 
pressed. For bosonic models, it is demonstrated that the causal region 
may be bent to the outside, covering distances exponential in time. 
Information can propagate from region A to B at arbitrary veloci- 
ties, via the quantum channel defined by the free time evolution of 
the chain. Note that this contrasts the situation in disordered spin 
chains, where the causal region is curved to the inside [9|]. 



rigorously, without resorting to numerical means. We do so 
by considering single excitation spaces and — in this context 
unusual — invoke ideas from convex optimization. 

There are several conceptual consequences of these results. 
It is now clear that any analysis of non-equilibrium processes 
in bosonic models must incorporate the possibility of far-away 
regions exchanging information on short time scales. In par- 
ticular, it seems likely that simulating short-term dynamics 
even of the low-energy sector of bosonic models is much 
harder than for spin chains (where Lieb-Robinson bounds are 
the basis for efficient algorithms la udO . Further, the results 
highlight the non-triviality of Lieb-Robinson bounds for spin 
chains with finite-dimensional constituents. 

More practically — while the models we present very 
strongly violate any bound on propagation speeds — they have 
reasonable physical properties. Related models with similar 
features could well be realized by tuning the parameters of 
suitable physical systems, for example in arrays of coupled 
cavities with polariton excitations [11]. This opens up the 
possibility of observing accelerating excitations experimen- 
tally and, potentially, of using bosonic chains as fast channels 
for quantum communication. 
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Local Hamiltonians and causality in spin chains. - A local 
Hamiltonian on n sites is of the form 

n 

H = Y,hj, (1) 
3=1 

where hj acts non-trivially only on a finite number of adja- 
cent sites. In what follows, we will restrict attention to the 
most relevant case of nearest-neighbor interactions. Quantum 
information transmission through spin chains with Hamiltoni- 
ans as above has been extensively studied in the literature. 

Before turning to bosonic models, let us first recall the 
precise situation for spin chains (d-level systems). The fact 
that there always exists a speed of sound — a maximal speed 
of information propagation — is the content of the following 
Lieb-Robinson bound [1]: If A and B are operators which act 
non-trivially only on some (distinct) regions of the chain, then 
there exist constants fi, C > and a velocity v > such that 

||[A(t),B]|| < C||A|| || B || e -M(di B t(A,B)-«|t|) ) (2 ) 

for all times t. Here, A(t) — e lHt Ae~ lHt is the time-evolved 
observable, ||.|| the operator norm, and dist(A, B) denotes the 
number of sites between the supports of A and B (see Fig. [T}. 
The above form may seem somewhat awkward at first sight. 
To get a more physical statement, one may verify that Eq. (O 
implies that any effect a perturbation A can have on a distant 
observable B is exponentially suppressed outside the causal 
cone defined by \t\ > dist(A, B)/v. In particular, any non- 
exponentially suppressed quantum communication using this 
spin chain can happen at most with velocity v JUl- Due to 
the intuitive nature of this statement the above bound is often 
taken for granted or even dismissed as being "trivial". 

Supersonic communication. - Roughly, we say that a model 
allows for "supersonic" communication, if its dynamics can 
carry information over distances m in time t{m) which scales 
sub-linearly in m. We will make this concept precise below. 

The setting is a chain of n bosonic systems with nearest- 
neighbor interactions and open boundary conditions. The in- 
teractions should be translationally invariant (hi = hj in Eq. 
(HJ), up to the obvious modifications at the boundary. We re- 
fer to the left sites 1, . . . , a as section A of the chain, whereas 
sites a + to, . . . ,n form part B. We assume that the system is 
initially in some factoring, translationally invariant pure state 
IV-') = IV'o)®"- A party in control of region A may now try 
to communicate with a party at B by either creating some ex- 
citations in her end of the chain, or else leaving the system 
untouched. More precisely, in the first case party A would ap- 
ply a unitary operator Ua to region A. At the receiving end, 
party B waits for some time t before probing whether a signal 
corresponding to some POVM element Ob is detected. The 
statistics are influenced by A's decision and given by 

P 1 =tx[0 B e- HH U A \^)(->P\Uy tH ] 
in case A has excited the chain and 

P =tv[O B e- itH \ 1 p){^\e itH ] 



in case A has not done so. The classical information capac- 
ity of the channel thus defined is a function of the signal 
strength 5 = \Po — Pi\. If 5 scales as 1/ poly (to), stan- 
dard protocols involving polynomially many channels used 
in parallel may be employed to, say, "transmit radio sig- 
nals through the quantum chain" with arbitrarily high fidelity. 
The last relevant quantity is the energy scale of the states in- 
volved, measured, e.g., by the variances Eq = (ip\H 2 \ip), 

Ef = (i/j\U^H 2 Ua\iP). Low values for Eq, E x imply that the 
states are largely contained in the low energy sector of H lfl3ll . 
Set e — max{i?o, -Ei}. Below, we define three increasingly 
strong ways in which bosonic models could potentially violate 
finite bounds on the maximum propagation speed for signals. 
We go on to establish the main result: even the strongest sce- 
nario can be realized by reasonable Hamiltonians. 

(i) Models which allow for arbitrarily fast transmission of 
information, using polynomial resources. More precisely, for 
every signal velocity m/t, there should be suitable encoding 
operations [/^(to) and observables Ob (to) such that the sig- 
nal strength 6(m) is of order 1 / poly (to). To obtain a reason- 
able protocol, the energy scale e(m) should grow only poly- 
nomially in in. While models of this type are interesting ob- 
jects of study, it may be argued that their existence would not 
be too surprising. Indeed, as energy and time take reciprocal 
roles in quantum mechanics, it is plausible that adding "more 
energy" to the system may lead to faster dynamics. This mo- 
tivates the next, more stringent, situation. 

(ii) Models for which the signal velocity scales faster than 
the inverse energy. In addition to the definitions above, we 
demand that mS(m) / (te(m)) — ► oo as m — * oo. For such 
models, the phenomenon cannot just be explained by the fact 
that unbounded Hamiltonians allow for signals with higher 
energies and thus faster dynamics. 

In scenario (i), (ii) above, information propagates at arbi- 
trarily high velocities — yet the distance covered is still linear 
in time (so the causal regions are cones with arbitrarily wide 
opening angles). The final situation is more demanding, re- 
quiring that excitations "speed up" as they propagate. 

(iii) Models allowing for accelerating signals. Here, we 
require that the signal strength 6, the energy scale e and, in 
fact, the encoding operation Ua do not depend on the distance 
to, while the time t should scale sub-linearly in to. In the next 
section, we discuss situations exhibiting behavior of type (iii) 
(and hence also of types (i,ii)) in a quite radical fashion. 

Models. - The type of models we subsequently allow for 
are governed by nearest-neighbor Hamiltonians of the form 

n— 1 n 
3=1 3=1 

with interaction term fj.j+i and on-site term gj . For fj,j+i = 
(aj + o>j)(Oj +1 + flj+i), and gj = fia^aj for fi > this is 
an instance of a harmonic chain. For an on-site interaction 

gj = iiOjCLj + Ua/jajfajCij — 1) and a hopping fjj+i = 

J(ajaj+i + h.c.) this gives rise to the Bose-Hubbard model. 

Let us spend some time to develop the physical intuition 
behind the constructions below. Consider a Hamiltonian with 
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harmonic hopping terms of the form fj.j+i — J(ajaj+i + 
/i.e.) and initial state — |0, . . . , 0} (using the Fock basis). 
Diagonalizing the operators fj,j+i, we see that the coupling 
strength between the sites grows as higher Fock layers become 
populated. Now consider an on-site interaction gj which does 
not preserve the Fock basis. If gj is e.g. a (low-order) poly- 
nomial in aj 7 Oj, then any application of gj will couple the 
state \ip(t)) to higher and higher Fock layers. For the sake of 
the argument, one may be tempted to model the spectrum of 
local reductions to undergo a random walk on the local Fock 
basis, starting at |0) and gradually spreading to higher lev- 
els. Therefore, for such Hamiltonians, the coupling strength 
between sites could be expected to grow as time proceeds. 
Terms coupling different Fock states are common in physical 
models (as long as the Fock basis does not represent massive 
particles). An obvious example is provided by the harmonic 
chain fjj+i — (a] + + Oj+i), whose Hamiltonian 

contains terms of the form aj<Xj+i. In particular, as Hamilto- 
nian dynamics by definition conserves energy, these terms do 
not inject energy into the system. 

It is unclear at this point whether this increase in coupling 
strength gives rise to an acceleration of the dynamics. Our 
objective below is to show that this can indeed happen, by 
constructing a model which exhibits an extreme violation of 
the usual causality bounds - allowing for exponentially accel- 
erating signals of constant strength. Other models compatible 
with the above intuition would naturally be expected to show 
similar accelerations, albeit not necessarily exponential ones. 
The model is constructed to violate the bounds in the strongest 
possible fashion while still being solvable. It serves as a proof 
of principle and as a worst case estimate for applications of 
Lieb-Robinson bounds for the simulation of dynamics. 

Specific "exchange interaction " model. - We will pay spe- 
cial attention to the following model, defined for bosons with 
spin 1 (so associated with the Hilbert space H = £ 2 (1R) ® 
(D 3 ). We define for site j the operators Aj-ki — \k, j | 
and Bj-kj = \k, I |. The Hamiltonian is specified by 



k,l=l 
oo 

g 3 = 2$^(ifc|*+l,T)(M | + /i.e.) + |0,I)(0,I |. 
k=i 

Note that fj.j+i may be looked at as a variant of the famil- 
iar exchange interaction. Clearly, H is a legitimate bosonic 
Hamiltonian with translationally invariant nearest-neighbor 
interactions. We will prove our claim three steps. 

1. Mapping to an excitation Hamiltonian. - To start with, 
|0, 1; . . . , 0, i) is an eigenstate of the Hamiltonian. If we now 
place a single excitation with spin f at the first site — so start 
with the initial state vector 1, |; 0, j; . . . ; 0, j) — we see that 
time evolution will only couple this to state vectors of the form 



110) 



|0,|;.. 
|o,4;-- 



I - i±i 1-0 I- 

,0,j;^2,T;0,i;. 



,0,1), if I odd, 
,0,|), if I even, 



with particles at site j = 1/2 + 3/2 and j = 1/2 + 1, respec- 
tively (for I = 0, . . . , 2n — 1). When considering only such 



excitations, we can hence pass to a new effective Hamiltonian 
E with Hilbert space K = <C 2n and initial condition 



E= J2 i(l + l)(ll* + l»«*H-||i»«* + l|l)- 

This Hamiltonian faithfully models the single excitation sec- 
tor in the above sense. Note that has now two roles: It 
both refers to a position in the original lattice, as well as the 
particle number in the original Hamiltonian. Also, to simplify 
notation later on, we will at this point pass to the half-open 
chain by setting n = oo 10]. The form of E makes it man- 
ifest that the excitation experiences stronger coupling coeffi- 
cients as it moves along the lattice. This alone is not sufficient 
to ensure an accelerated, directed movement: convincingly, 
the ramping up in the coupling coefficient could cause the 
excitation to scatter back. Maybe surprisingly, we will find 
travelling solutions below. 

2. Moments. - In this single excitation sector, time evolu- 



tion corresponds to p(t) 



,-iiJS 



p(0)e itE , for states on JC. It 



proves expedient to introduce the operators 



1=0 1=1 

Note that X corresponds to a discrete position operator, mea- 
suring twice the distance of the original model. There is a lot 
of structure in this model: the commutation relations between 
these operators form a closed algebra. Indeed, 

i[E,X] = P, i[E,P] = AX - 21, 

the algebraic completion of which being that of the Lie- 
group su(2). Using the familiar Baker-Hausdorff formula, the 
Heisenberg picture time evolution of X under E is given by 



X(t) 



e itE Xe~ itE 



X + t[iE,X} + — [iE, [iE,X]]- 



One obtains a closed-form expression for X(t) by exploiting 
the relations in the algebra A — {E, X, P, 1} to iteratively 
solve these nested commutators. Explicitly 



oo 

X(t) = 1 + ^- 



-,22-142-1 

>i<»--W H "- n 'ti <*>< 

X - P^sinh(2t)+(X - - 1) (cosh(2t) - 1), 



For the time evolution of X starting from the single excitation 
||0)), we hence find 

((0|Jf(*)|0)) = i(l + cosh(2t)). 

Thus, the expectation value of X is increasing exponen- 
tially in t. This fact alone, however, is not enough to show 
that we have signaling: It could be that the excitation de- 
velops a long asymptotic tail that leads to large first mo- 
ments, but carries a small weight. Thus, further information 
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is needed. It turns out that knowledge of the second mo- 
ments ( (0 1 |X (i) 2 1 1 0) } - ( (0 1 1 X (f ) 1 1 0} } 2 is sufficient to prove 
signaling using the convex optimization ideas below. An 
analogous — if more tedious — calculation arrives at 

((0\\X(t) 2 \\0)) = cosh 2 (t)cosh(2t). 

3. Hitting time from a convex optimization problem. - For 
some site m > 1 define the hitting operator 

oo 

I=2ra-1 

So Pi = tr[p(t)T] is the signal a distant observer m sites away 
from the origin may receive, compared to Po = 0. Hence, 
the set A = {1} is the single first site, whereas B — {m + 
1, . . . , oo} is the natural right part of the chain. Let us set 
M = 2m — 2. We will now bound this expectation value by 
analytically solving a convex optimization problem: 

oo 

minimize pi , 

l=M+l 

oo oo 

subject to YVZ + l)pi = a(t), VVZ + l) 2 Xi = b(t), 

1=0 1=0 

YaZoPi = 1> an( ^ Pi — ^ f° r a H where a(t) = (1 + 
cosh(2t))/2 and 6(f) = cosh 2 (t) cosh(2f). So we minimize 
the signal, given first and second moments, as a worst case 
analysis. This is an (infinite-dimensional) linear program, and 
hence an instance of a convex optimization problem. We can 
readily get a bound to the optimal solution by identifying a 
suitable solution to the Lagrange dual problem 11511 . which is 

maximize —d T y, 
subject to F T y > —c, 

where F ld = j, F 2 j = j 2 , F 3J = 1 for j = 1,2,..., and 
d(t) = (a(t),b(t),l) T . Also, c = (0, . . . , 0, 1, 1, . . . ) T , as 



a vector starting with M zeros. The latter is an optimization 
problem over y. Any solution we can identify of the dual will 
give a lower bound to the primal. It can be shown explicitly 
that the subsequent vector y* = (-2(1+M)/M 3 , 1/M 4 , (1+ 
M) 2 /M 2 — 1) T is always a feasible solution of the dual prob- 
lem. Taking the time t* — log(M) gives 

c T x* > -d T y* = 2a(log(M))(l + M)/M 3 

- 6(log(M))/M 4 - (1 + M) 2 /M 2 + 1 =: g(M). 

For this function we have that limM->oo g{M) = 3/8, and 
g{M) > 1/5 for M > 9. The solution of the dual will hence 
give rise to a lower bound of our primal problem. Hence, af- 
ter a time t* logarithmic in m, a signal of constant strength 
tr[Tp(f*)] > 1/5 will have reached party Bl This means, of 
course — within the validity of the model — that we can in prin- 
ciple signal at any speed over arbitrary distances: The signal 
will not even decay, and the Holevo-x and the classical infor- 
mation capacity of the associated quantum channel are indeed 
constant. One can communicate with an exponentially accel- 
erating signal of type (iii) in the above classification. 

Summary and Outlook. - In this work, we have shown that 
in systems of locally interacting bosons, excitations may ac- 
celerate and carry information faster than any finite speed of 
sound. Many results on the simulatability of dynamics and 
ground-state properties of spin systems (in particular using 
DMRG techniques) rely on Lieb-Robinson bounds 11611 . Ex- 
ceeding the usual difficulties associated with reasoning about 
infinite-dimensional systems, our examples imply further dif- 
ficulties any generalization of such results to bosonic systems 
has to deal with. This remains true even if one restricts atten- 
tion to low-energy sectors. Interestingly, while analytical and 
numerical treatments of such models face formidable difficul- 
ties, they may be simulated using analogue physical systems, 
such as coupled cavity arrays 111 ill where massless polaritonic 
excitations are not subject to particle number conservation. 
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I. APPENDICES 

For clarity of the argument, we present a few aspects of 
the above presentation in some more detail. Note that this 
material is identical with the above one. 



This is obviously an eigenstate of the Hamiltonian. Now we 
excite site 1, forming set A = {1}, by setting it to 1, |), while 
keeping the rest of the chain unaltered. We then see that the 
initial state vector after the excitation 

couples to 

|o,j)|U)|o,i)^"- 2 \ 

and this to the former and to 

|o,i)|2,T)|o,j)^"- 2 ). 

Now the pattern is clear: All excitations are contained in 
span {|0,i)® fc |fc,i)|0,i)^"- fc - 1 \ 

|oa)®1* + i,T)|o,J,)® (n ~* -1) }, 

and each basis vector — except at the ends of the chain — 
couples to two further basis vectors. It is now a straightfor- 
ward exercise to verify that in this subspace of excitations, 
one arrives at the above given effective Hamiltonian. 



Appendix A: The spin Hamiltonian 



Appendix B: Convex optimization problem 



We consider a bosonic model with a spin-1 internal de- 
gree of freedom. The local Hilbert space is hence spanned 
by {|fc, t), \k, i) : k £ N }. We will be able to keep track of 
local excitations in the following model. We take the Hamil- 
tonian 

n , oc 

h = 5i;(5^(2*-i)(i|fc,i><i,Ti®ii,4.)<fc,ii+fc.c.) 

j=i V,;=i 
+ 9 j 

where the on-site interation is taken to be 

oo 

9j = 2^(ifc|fc + 1, T)(fc, 4 1 + h.c.) + |0, |)(0, 1 1. 
fc=l 

The hopping term is a variant of an exchange interaction be- 
tween neighboring sites (except from the spin, this would ex- 
actly be a standard exchange interaction). This Hamiltonian 
can also be written as 



n / oc 



H = E( E ( 2fc - 1 )( iA h,k B *+w> + l, - r - ] ■ ! <* 

3=1 ^k,l=l 

with the hopping operators at site j being defined as 

^• J k,j = IM>(i»H B j]ktl = \k,i)(i,l\, 

The initial state vector is taken to be 

iv>> = io,ir n . 



The matrix F in the convex optimization problem defining 
our hitting time problem is given by 



F = 



12 3 4 
1 4 9 16 
1111 



The first row captures the first moments, the second the 
second moments, whereas normalization of the probability 
distribution is enforced by the third one — together with the 
fact that all entries of x are positive. The vector c = 



(0, 



,0,1,...,!), as a vector starting with M zeros. The 



primal problem can hence be written as 

minimize c T x, 
subject to Fa; = d, 

x>0, 

where d(t) = (a(t),b(t), 1) T . This is a linear program in 
standard form, with a matrix inequality constraint, and posi- 
tive elements of the objective vector x. In this form it is easiest 
to identify the dual problem. It is given by 

maximize —d T y, 
subject to F T y > —c, 

where now y is not constrained to be positive. The constraints 
on the vector y can equally be written as 



jyi + 3 2 V2 + y-i + 8j>M > 



(3) 



for all j, where Sj>M takes the value 1 if j > M and is zero 
otherwise. The question is: Can we find for each M a choice 
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of (j/i, 2/2, 2/3) and a time £ such that we can meaningfully find 
a positive solution to the dual problem? 

We can take any solution to the dual problem. We choose 

Vl = -2(1 + M)/M 3 , 
V2 = 1/M 4 , 
2/3 = {l + Mf/M 2 



1, 



and will verify that it is indeed a solution with the appropriate 
properties. The intuition behind this construction is as fol- 
lows: Consider the quadratic function / : [0, 00) — > 00 as 

f(x) = xy 1 + x 2 y 2 + 2/3. 

This is identical with the left hand side of Eq. ||3}, up to being 
defined on [0, 00). We require that / takes the value for the 
first time exactly at x = M, and that f(x) > — 1 for all x. 
Taking 2/2 = 1 /Af 4 then gives rise to the above construction, 
satisfying Eq. (0 for all M and all j. We have jyi + j 2 y 2 + 
2/3 = exactly at j = M, which is a desirable feature to get 
a meaningful bound to the original problem at hand. This is 
hence always a feasible solution to the dual problem. We will 
therefore obtain a lower bound to the optimal solution of the 
primal problem c T x* as 



T * 
C X 



> -d T y 

= 2a{t){\ + M)/M 3 - b(t)/M 4 
- (l+Mf/M 2 + 1. 



Choosing the time t* 
we find 



log(M), logarithmic in the distance, 



T * 
C X 



> 2a(log(M))(l- 
- (1 + M) 2 /M 2 



M)/M 3 - &(log(Af))/Af 4 
- 1 =: 9(M). 



Now it is not difficult to see that indeed, in the limit of large 
M we have 



M 



lim g(M) = 3/8, 



and that g(M) > 1/5 for all M > 9. This means that for 
each large m — and hence each large M — we can arrive at a 
constant signal for an appropriate constant signal. 



Appendix C: Bose-Hubbard-type models 

In a sketchier fashion, we will now briefly address the ques- 
tion whether even Bose-Hubbard-type models 



particle at site A — {1} compared to \N)® n — will lead to 
a signal at site B = {n} that is at most polynomially sup- 
pressed in n. The signal increases linearly in N. Hence, hop- 
ping models can also display violations of any finite bound 
on propagation speeds, at least in the sense (i) of the above 
classification. 

For each chain length n, we will consider a relevant n- 
dimensional subspace of the original Hilbert space 



C 



N 



span{|ji 



,j n )\3k:j k = N + l,ji = NVl^k}, 



for a suitable integer TV, and denote with = H\£n 
the orthogonal complement of the infinite-dimensional Hilbert 
space. The initial state vector of the chain of length n is then 
taken to be 



W) = \N + 1)\N) 



®(n-l) 



e C 



N- 



This means that we fill up all number states with bosons at 
each site up to N, and add an additional particle at site A = 
{1} as a local excitation. We will now see that we can then 
approximate the dynamics of the chain within a spin chain 
restricted to a single excitation Cn arbitrarily well. Define for 
site j the operator Aj = \N){N + 1|, then the above chain is 
approximated in its dynamics by the dynamics of 
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V=(N + l)Y / (A]A J+1 +h.c. 



This dynamics is specifically simple, and we can see that in 
a time t* that is linear in N we will receive a signal at site 
B = {n} that is at most polynomially suppressed in n. This 
clearly means that this signal will be received at an arbitrarily 
velocity, if ./V is sufficiently large. To see that this mapping 
can be done to arbitrary approximation, we need the subse- 
quent observation: 



Lemma 1 For any Hermitian matrix M, partioned as 

M 



A B 

B^ C 



• have that for all t 6 [0, 00), 

lim sup \\Oc itK Oi -e itA \\ =0 



where 

M x = I M 




P 



P>xl\, 0=[t 0] 



H = ^2 ( ( a ] a j+i + h.c.) + h(a,jO,j) 



Here, ||.|| denotes the operator matrix or vector norm. To 
prove this, note that 



with some function h : R + — > R + can display such a behav- 
ior. Indeed, it can be seen that if h increases in n sufficiently 
strongly, leading to highly interacting particles, then there ex- 
ists for each chain length n an N (polynomial in n) such that 
the initial state |?/>) = \N + l)\N)®( n - 1 '>— so an additional 



d 



A B 
fit c + P 



= inftr 



A B 

St C + P 



where the infimum is taken over Q that are projectors of rank 
d, when A is a d x c?-matrix. This infimum exists for any 



7 



M E M x , call it Qm- Since \\A\\ = c is constant, we have 
that 



lim sup 



1 




= 0. 



Hence, since \e ltx i ^ | = 1, the assertion follows. 

Subsequently, the submatrix A will be identified with H 
restricted to Cn, whereas C is the restriction to the orthogonal 
complement Cj^. Using this observation, we see that for each 



e > we can find a function h : 1R + — > R + such that the 
following property is satisfied: For each chain length n there 
exists an N such that for the above hitting time t* 



-it'Hi 



-WV 



< e 



for all t e [0, t*]. But this means that we can signal with this 
chain, in the sense of (iii) in the above classification, as the 
class of restricted models has this property. 



